Journal of Statistical Physics, Vol. 86, Nos. 5/6, 1997 

Corner Exponents in the 
Two- Dimensional Potts Model 
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The critical behavior at a corner in two-dimensional Ising and three-state 
Potts models is studied numerically on the square lattice using transfer op- 
erator techniques. The local critical exponents for the magnetization and the 
energy density for various opening angles are deduced from finite-size scal- 
ing results at the critical point for isotropic or anisotropic couplings. The 
scaling dimensions compare quite well with the values expected from confor- 
mal invariance, provided the opening angle is replaced by an effective one in 
anisotropic systems. 
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1. INTRODUCTION 

The surface shape of a system may have some influence on its local critical 
behavior at a second order phase transition. This was shown by Cardy^ 
for a magnetic system with O(N) symmetry within mean-field theory and in 
d = 4 — e dimensions. The local magnetic exponent at a wedge was found to 
vary continuously with the opening angle 9. 

Marginal behavior was also obtained at a corner in the two-dimensional 
Ising model using the star-triangle recursion relation on the triangular lattice 
and the corner-to-corner spin correlation function on the square lattice to cal- 
culate the corner magnetization. ( 2,s ) The same problem was later studied on 
the square lattice using row-to-row^ 4 ) or corner® transfer matrix techniques. 
An expression for the 90° corner magnetization was conjectured in ref. 4. 
Recently, Abraham and Latremoliere^ 6, 7 ) obtained the edge magnetization 
analytically, as a function of the distance from the corner, thus confirming 
the conjecture of Kaiser and Peschel. 
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Other systems were also considered, such as the planar Heisenberg anti- 
ferromagnet^ 8 '' and the self- avoiding walk confined into a wedge in two and 
three dimensions. ( 9 ~ 13 ) 

In two dimensions a varying corner exponent x c (8) immediately follows 
from the conformal mapping w = z 6 '* which transforms the half-plane into 
a wedge with opening angle 9S 14 ' 2 ^ The decay of the critical corner-to-bulk 
correlation functions then gives 

x%9) = ^x s (1.1) 

where x s is the corresponding surface exponent. This result is valid for 
isotropic systems only. When the couplings are anisotropic, lengths have to 
be rescaled in order to restore isotropy 1 - 1 ' 2 ** and (|1.1| ) still applies with an 
effective value of the opening angle. 

The marginal local critical behavior may be understood by considering a 
system with a "parabolic" shape, y = ±Cx a S 15 ' 1& > Under a length rescaling 
by a factor b, C transforms into C = b a ~ 1 C. Thus the flow is towards the flat 
surface geometry when a > 1 and towards the half- line when a < 1. When 
a = 1, i.e., for the corner geometry, the surface shape is scale invariant and 
C = tan(#/2) is the marginal variable. 

Up to now, the conformal result ( |1 . 1| ) has been checked for the local mag- 
netic exponent x c m {9) with opening angles corresponding to simple fractions 
of 7T. In the present work, we extend these results by studying the local criti- 
cal behavior of both the energy and the magnetization in the two-dimensional 
g-state Potts model with q = 2,3. Various rational values of tan(#/2) and 
tan# are considered. In Section |2| we use a local operator formalism for the 
construction of the row-to-row transfer matrices which are needed for the 
corner geometry on a square lattice. The corner exponents are deduced from 
a finite-size scaling analysis of the data for isotropic and anisotropic systems 
in Section [3]. 

2. POTTS MODEL IN THE CORNER GEOMETRY 

We consider the zero-field Potts model with Hamiltonian 
- /3H = Yll K i (<7<^(i,jM*+i, j) - 1) + K 2 (q$*(i,j),a(i,j+i) ~ 1)] (2.1) 
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Fig. 1. Corner geometry and transfer operator. 



where 5 is the Kronecker delta function and the Potts variables a(i, j) take 
the values 0, 1, • • • , q — 1. When q = 2, the Ising model is recovered. The 
sum runs over the bonds of a square lattice in the corner geometry as shown 
in Fig. 1. The couplings are assumed to be anisotropic with values K\ in the 
vertical direction and K 2 in the horizontal one. 

Let us define the row-to-row transfer matrix elements as 



T m ,m—1 v7t|^m|^ 1) 

N(m) 

exp 



N(m)-1 

j=n(m) j=n(m) 



(2.2) 



where Oj (a'j) denote the Potts variables on the mth [(m — l)th] row, respec- 
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tively. \m) is a state vector corresponding to a configuration of the Potts 
variables from n(m) to N(m) along the mth row and T m is the transfer 
operator acting on this state vector. 

The transfer operator can be written as the product T m = V2(m)Vi(m) 
with ( 17 ) 



VAm) 



Vo(m) 



cxp 



cxp 



N(m) <j-l 

k{ E E(^) 

j=n(m) p=l 
JV(m)-l 5-1 

K2 E E(c]^ + i) p 

j=n(m) p=l 



(2.3) 



K\ is a dual coupling such that x\ x\ = q with Xi = exp(gKj) — 1. The 
operator Cj is diagonal on the basis of the Potts states \aj) whereas Rj acts 
clS cL ladder operator on the same basis, so that 



Cj\<7j) = exp (i — Uj^j \(7j) , (Rj) p \<jj) = \<jj + p) (mod q) 



(2.4) 



Introducing local transfer operators 



<?-i 



h(j) =xr + E(^) P , t 2 (j) = 1 + E(^j^ + i) P 

p=0 p=0 



9-1 



(2.5) 



we have that, up to an unimportant constant factor, the row-to-row transfer 
operator appearing in ( |2.2| ) can be rewritten as 



N(m)-1 N(m) 

T m = II HO) II Hi) 

j=n(m) j=n(m) 



(2.6) 



This transfer operator may be viewed as a discrete time evolution operator 
acting on the Potts state vector at time m — 1 to make it evolve to its state at 
time m. Due to the corner geometry, this operator is time- dependent. Given 
an initial state \i) at time 0, corresponding to the top edge in Fig. 1, this 
state will evolve to 

\m) =T m T m _ 1 ---T 2 T 1 \i) (2.7) 

at time m. In the sequel we shall use either free or fixed boundary conditions 
on the top edge. Fixed boundary conditions correspond to an initial state 
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I A;) = \u1U2 ■ ■ '&n) whereas for free boundary conditions one has to take a 
superposition |/) = q~ N l 2 YX=\ 1^) °f the Q N configurations of the N Potts 
variables with equal amplitudes. One may notice in Fig. 1 that a site variable 
Oj with j < n(m) or j > N(m) remains frozen at time m and later. 

The corner magnetization m c is calculated with the Potts variable of, 
where t is the index of the column corresponding to the tip (see Fig. 1). It 
is given by^ 18 ^ 

- 1 , 1 (f\C't\M) 



where the average of the Kronecker delta follows from (|2.4j) . The state |/) cor- 
responds to free boundary conditions on the sides of the wedge whereas \M) 
is defined as in (|2.7| ) with \i) = 1 00 - - - 0) . Thus the top edge is fixed with all 
the sites in the state 0. At the critical point this symmetry-breaking bound- 
ary condition ensures that the tip magnetization remains non-vanishing on 
a finite system. 

The energy density associated with a bond can be defined, up to a con- 
stant factor, as the average of the Kronecker delta. When the bond is hor- 
izontal, the system evolves as above up to time M — 2, where the value of 
S at -i,(Tt i s taken (see Fig. 1), and the corner energy density is given by 



h ls \ l q ^{f\T M T M -i{CUC t y\M-2) 



e 



c 



With a vertical bond one obtains nonvanishing contributions to the average 
of the Kronecker delta when there is no flip along this bond, so that ti(t) 
in ( p.5|) contributes a factor 1 + x\. With the geometry of Fig. 1 the energy 
density can be written as 

^MMiN ^-'^y^ (2.10) 

These expressions are easily generalized for other shapes. 

3. NUMERICAL RESULTS 

The Potts exponents for the corner magnetization and the corner energy 
density have been obtained through finite-size scaling at the critical point for 
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Fig. 2. The different system shapes used to calculate the corner magnetization and 
energy density. The Potts states at the top eddge (bold line) are either fixed or free. 



different opening angles, using the shapes sketched in Fig. 2. In order to limit 
the size N of the system, we use only unit steps (p = 1) in the horizontal 
direction and change the angle 9 by varying I (see Fig. 1). The maximum 
sizes studied are N = 19 for q = 2 (Ising model) and iV = 11 for q = 3, so 
that an extrapolation of the finite-size estimates of the exponents is needed. 
Isotropic and anisotropic sytems have been considered. 

When the system is isotropic the critical point, corresponding to x±x 2 = q, 
is located at K c = (1/q) ln(l + y/q). The critical corner magnetization decays 
as 

mf\N) = A m (9) N~ x ™ 

with the system size N. 

The corner exponent x^iO) deduced from sequence extrapolations of the 
finite-size estimates using the BST algorithm* 19 ) is given in Table I. The 
numerical results are in quite good agreement with the values expected from 
conformal invariance in Eq. fll.lj) with the surface magnetic exponents at the 
ordinary surface transition x s m = 1/2 for q = 2 and x s m = 2/3 for q = 3.( 14 > 16 ) 
The critical corner energy density 

ef 1 (TV) = ef *(oo) + A e {6) N~ x * 
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Table I. Scaling Dimension of the Corner Magnetization for the q-State Potts Model 
as a Function of the Opening Angle." 



(deg) 






(deg) 






q = 2 


Expected 


q = 3 


Expected 


360(d) 


0.25(1) 


0.25 


180(e) 


0.667(1) 


0.6667 


270(c) 


0.332(1) 


0.333333 


90(b) 


1.331(3) 


1.3333 


225(c) 


0.400(1) 


0.4 


53.13(b) 


2.243(3) 


2.2586 


180(e) 


0.5000(1) 


0.5 


45(a) 


2.66(1) 


2.6667 


90(b) 


0.999994(6) 


1 


26.56(a) 


4.516(3) 


4.5172 


53.13(b) 


1.693961(7) 


1.693955 


18.43(a) 


6.52(2) 


6.5094 


45(a) 


2.00000(2) 


2 


14.04(a) 


8.53(4) 


8.5493 


36.87(b) 


2.44098(6) 


2.441016 


11.31(a) 


10.58(4) 


10.6101 


28.07(b) 


3.2057(3) 


3.205986 


9.46(a) 


12.66(4) 


12.6819 


22.62(b) 


3.979(4) 


3.978804 








11.42(b) 


7.8796(6) 


7.880092 









"The numbers in parentheses give the estimated uncertainty in the last digit. The letters 

refer to the system shapes in Fig. 2. 



contains a regular part which depends on the opening angle and the bond 
orientation. In order to obtain the corner exponent the regular contribution 
has to be substracted. This can be done by calculating e^ nt (iV) for either 
fixed or free boundary conditions at the top edge. The asymptotic value is 
the same in both cases, but the amplitudes of the finite-size corrections are 
different. Taking the difference, one obtains 

Ae c (iV) = [A^ xcd (9) - A^ cc (d)} N~ x ° 

Proceeding in this way, one avoids a systematic error, linked with the estima- 
tion of the regular part, and the finite-size corrections are amplified because 
the signs of the amplitudes are different for the two boundary conditions. 

The extrapolated values of the finite-size estimates for the corner expo- 
nent x c e {9) are given in Table II. They do not depend on the bond orientation, 
as expected, and are identical for q = 2 and q = 3. This is consistent with 
the conformal result, since, at the ordinary transition, the surface exponent 
x s e which enters in is equal to 2, the dimension of the system, quite 
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Table II. As in Table I, for the Corner Energy Density. 



(deg) 




x e 




q = 2 


q=3 


Expected 


180(e) 


2.001(1) 


2.01(3) 


2 


90(b) 


3.999(1) 


4.00(4) 


4 


45(a) 


7.993(7) 


8.0(1) 


8 


26.56(a) 


13.52(2) 


13.5(1) 


13.5516 


18.43(a) 


19.6(4) 


19.7(4) 


19.5281 



generally. ^ 20 ' 16 - ) The numerical values are in good agreement with the results 
expected from conformal invariance, too. 

The conformal expression for the corner exponent ( |1.1[ ) only applies to 
systems where the correlations are isotropic. In the case of an anisotropic 
system, with a coupling constant ratio r = Ki/K 2 ., the correlation lengths 
are different and take the values £i and £2 in the vertical and horizontal 
directions, respectively. Isotropy can be restored by changing the lattice pa- 
rameters di and 02 in such a way that, in the rescaled system, the correlation 
lengths become the same, i.e., £iai = ^zazS 1 ' 2 ' As a consequence, one obtains 
an effective opening angle which, in the geometry of Fig. 2a, is given by 

tan 9 cS = ( tan 9 , £ = — = — ^ (3.1) 

where ( is the anisotropy factor. 

For the Potts model, the following form of the anisotropy factor at the 
critical point has been conjectured: ( 21 ** 

> , ttm sinw x lc ^/q 

C = tan 2A' imTA^) = 7r^' 2cosA = v^ (3-2) 

When q = 2, it reduces to the known exact result ( = cosh2i^i c / cosh2i^ 2c 
for the Ising model. ^ 

The corner exponents x c m and x c e obtained on anisotropic systems with an 
opening angle 9 = 45° in Fig. 2a and different anisotropy ratio r are given in 
Tables III and IV. They are in good agreement with the conformal result (|1 . 1| ) 
when the opening angle is replaced by its effective value given by (|3.1|). 
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Table III. Corner Exponents for the Ising Model on an Anisotropic Lattice as a 
Function of the Anisotropy Ratio r. a 



r 


numerical 


Expected 


Numerical 


Expected 


0.200 


4.5737(8) 


4.57374 


18.2(6) 


18.295 


0.300 


3.595(2) 


3.59518 


14.36(8) 


14.381 


0.400 


3.05(2) 


3.06772 


12.24(8) 


12.271 


0.500 


2.7326(1) 


2.73264 


10.92(5) 


10.931 


0.600 


2.4999(5) 


2.49960 


9.96(5) 


9.998 


10.00 


1.16(2) 


1.16200 


4.6(1) 


4.648 



11 We used the shape of Fig. 2a with an opening angle 9 = 45° . The numerical values are 
compared to the conformal result with a rescaled angle, as explained in the text. 



Table IV. As in Table III, for the Three-State Potts Model. 



rn 



r 


Numerical 


Expected 


Numerical 


Expected 


0.109 


9.31(3) 


9.30667 


27(1) 


27.920 


0.166 


6.98(4) 


7.00928 


21.2(9) 


21.028 


0.259 


5.30(3) 


5.31052 


15.9(2) 


15.931 


0.414 


4.05(4) 


4.05604 


12.1(1) 


12.168 


0.692 


3.137(6) 


3.13157 


9.3(1) 


9.365 


1.241 


2.45(1) 


2.45244 


7.3(1) 


7.357 


2.561 


1.955(6) 


1.95572 


5.8(2) 


5.867 


7.530 


1.5(1) ' 


1.59445 


4.6(2) 


4.783 
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4. CONCLUSION 



The scaling dimensions of the tip magnetization and energy density at 
a corner have been calculated for the Ising and three-state Potts models in 
two dimensions as functions of the opening angle 9. Using transfer operator 
techniques and finite-size scaling at the critical point, we have tested the 
result of conformal invariance, predicting a simple relation between corner 
exponents, surface exponents, and opening angle, for a broad spectrum of 
6 values. In the case of the Ising model, the magnetization results extend 
previous studies where simple fractions of tt were considered, whereas the 
energy results are new. Isotropic and anisotropic systems have been treated. 
In the later case, the conformal result still applies, provided the opening 
angle is properly rescaled in order to restore isotropy. 

As a possible extension of this work, one may mention the possiblity 
to consider the same problem for noninteger values of q, using the relation 
of the Potts model with either Withney polynomials'- 22 -' or ice-rule vertex 
models. ( 23 < 24 ) 
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